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I. INTRODUCTION 



It is by now a well-appreciated fact that non-Abelian gauge theories display many inter- 
esting and beautiful properties which are, no doubt, indicative of their deeply geometrical 
nature and symmetry. Already at the level of perturbation theory, many of the scattering 
amplitudes, in Quantum Chromodynamics (QCD) for example, have a simple and elegant 
form, although large numbers of Feynman diagrams have to be summed up to arrive at 
these PJ. This feature was originally understood in terms of calculational techniques based 
on string theory as well as recursion rules B- [ffl. (A field theoretic understanding of the 
string-based techniques has also emerged ||.) There are also indications that there are in- 
tegrable field theories hidden in non-Abelian gauge theories which describe some aspects of 
these theories, such as scattering for certain specific choices of helicities as well as scattering 
in the Regge regime 0- ||. Given these features, it is clear that the theoretical exploration 
of the structure of even the perturbative scattering matrix can be quite useful. 

Within a field-theoretic approach, recursion rules for scattering amplitudes have been 
very useful in understanding color factorization and many other properties. These rules 
were originally derived for tree-level amplitudes using diagrammatic analyses by Berends 
and Giele and have later been extended upto the one-loop level []TT|- [|12|]. Elaboration 
and extension of this technique, we feel, can be very fruitful. In this letter, we present a 
derivation of the recursion rules within a functional formalism without reliance on diagram- 
matic analysis and in a way valid for any field theory. From the basic equations, recursion 
rules valid upto arbitrary loop order can be obtained, although at the expense of increas- 
ing algebraic complexity. We explicitly display the tree-level and one-loop recursion rules. 
Renormalization constants, which one must consider beyond the tree-level are also easily 
incorporated in a functional derivation. 
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II. THE S-MATRIX FUNCTIONAL 



We start by considering a scalar field theory with action of the form 

S[4>] =S [<f>] + SU<j>], (2.1) 

where So[<f>] is the free action of the form So[(p] = J ^(pKcp and Si nt is the interaction part 
of the action. Specifically, K can be taken to be of the form K = —Z 2 (d 2 + m 2 ). The 
functional for the S-matrix elements can be written as |H| 



Ffy] = e F J [d(j>] e «*>M+iSi»*(*+*0 ; (2.2) 
where <p(x) can be expanded as 

= J2 a k u k(x) + a* k u* k (x) . (2.3) 
k 

u k (x),u k (x) are the one-particle wave functions and are solutions of the free field equation 
Kip = 0. F is given by Ylk a k a k an d is introduced in ( |2.2p to give the matrix elements 
for subscattering processes where some of the particles fly by unscattered. The matrix 
element for a process where particles of momenta k%, k% . . . scatter to particles of momenta 
Pi,P2- ■ ■ Pm is given by 

5 S 5 5 5 5 . . 

- — — F[tp] . (2.4) 

5a kl 5a k2 da kN da* pi da* 2 5a* M J^ =Q 



S, 



fel,fc2...fciV— »P1P2— PM 



Since F is not particularly relevant for our calculations below we can drop it in what follows. 
Further we consider (p(x) to be an arbitrary function rather than a solution to Kip = 0. 
Eventually, in obtaining the S-matrix elements we can choose it to be a solution to Kip = 0. 
We thus define 



\d(t>] eiSW-if^ (2.5) 



where we added a term exp(— i/2 J ipKip) for simplifications in what follows; this of course 
does not contribute when Kip = 0. In terms of this J 7 , the S-matrix elements are given by 



Ski , k 2 ■■ .k N ~>p 1 p 2 . . -pm 

where 



a kl . . . a kN a* pi . . . a* pM ^[ip}\ ^ . (2.6) 



at = /^4 (2 - 7) 
From the second of equations ( |2.5|) , we see that T is the generating functional of the con- 
nected Green functions with source Kip and (|2.6|) represents the well-known LSZ reduction 
formula. 

Consider now Differentiating with respect to ip, we find 

iK- 1 ^ = I \d<j>] ^sm-ifvK^ (2 8) 

din J 



Effectively, 0's inside the integral behave as iK j-. From Q2.5| ) we also have the equation 
of motion 

~5S 



-Kip 



J,S[<p]-i J <pK(f> _ (2.9) 



which can be written as 



i-r- + P($)F ~ (K<P)F = ■ (2-10) 
dip 

where p{<f>) = (^) and = zi^ 1 ^. 

J- itself generates connected as well as disconnected scattering processes. If we write 
T[ip\ = e lC ^\ where C[ip\ describes connected processes only, (p.lOj) reduces to 

- — C = Kip- T- x pQ>)T 
dip 

= Kip - P(-K~ l6 ^ + }) . (2.11) 

This can be regarded as a nonlinear functional differential equation for the S-matrix gen- 
erating functional and can be used for deriving systematic recursion rules for scattering 
amplitudes. We shall do this for Yang-Mills theory in the next section. 

Although it is not crucial to the discussion of recursion rules, one may also work with 
the quantum effective action or the generating functional for one-particle irreducible vertices 
r[$] defined by 



r[$] = c[<p] + J {Kip)®, 



with the connecting relations 



dip 



Then, ( 2.1 1Q becomes 



K§ = Kip - p($ + 0)1 



(2.12) 



(2.13) 



(2.14) 



or, in terms of 



5T 
5$ 



AT$ + p($ + 0)l 
fSS\ 



(2.15) 



=$+iA r - 1 ^- 



The right hand side involves the derivative of <& with respect to ip which is the two-point 
correlator G given by 



G^y)= fiff" 1 ^! S(y) 



(2.16) 



It satisfies the basic Schwinger-Dyson equation for the theory, viz. 

5 2 T 



G(y,z)=i5^\x-z). (2.17) 

Equation fl2.15|) supplemented by (|2.17|) can thus be regarded as a nonlinear equation for 
r[$]. From the connecting relations ([2.1 3|) we find that |^ = for Kip = 0, as is appropriate 
for S-matrix elements. In this case, from the definition of T, we have 



c[<p] = r[*\u s 



o • 



The S-matrix is then given by 



(2.18) 



(2.19) 



This relation gives a nonperturbative definition of the S-matrix. The free data in the solu- 
tions to |J = are the quantities on which T depends. (Perturbatively, the free data are 
the amplitudes a k , a* k in the solution for ip.) The fact that the S-matrix can be obtained as 
the exponential of (i times ) the action evaluated on solutions of the equations of motion is 
rather well known ||10|| , [ |i~3[ . 



III. THE S-MATRIX FUNCTIONAL FOR YANG-MILLS THEORY 

We start with the gauge-fixed Lagrangian £ of an S"?7(iV)-Yang-Mills theory given by 

C = -\F^ a - 1(0 • A) 2 - c(-d ■ D)c 

= ^Tr F^F^ + Tr (d ■ A) 2 + 2Tr c(-d ■ D)c, (3.1) 

where 

Ffj, u = F* u T a = d^Ay - d v Ap + g[A„ A v ] , 

DpC = d li c + g[Ap,c\ t (3.2) 

with [T a ,T b ] = f abc T c and Tr T a T b = — |5 a &. T a are matrices in the fundamental repre- 
sentation of SU(N). The free part and the interaction parts of the action are respectively 
identified as 



C = Tr A^-d^A" + 2Tr c{-d 



2 V 



CSt = 2gTr d,A u [A^ A"} + ^Tr [A„ A,} 2 (3.3) 
42 = 2^Tr^c[^,c]. (3.4) 

Following the discussion of Sec. 2, the functional T which gives the S-matrix can be con- 
structed as 

T\a\ = J e iS ~^ a ^~ d2 ^ a [dA}[dc][dc} , (3.5) 

The transition amplitude for N gluons of momenta ki, polarizations and colors labelled 
a\, . . . , qn to go to M gluons of momenta pj, polarizations e$ and color labels aw+i, • • • , % 



is 



From ( |2.10|) T satisfies the equation 



. (3.6) 

a=0 
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6a° ^ J 6a° 

H 1 H 1 



d 2 al? + J;(A)\ F+l 2gTr d,c[T a , c]e iS , (3.7) 

fin. " 



where 



ro(l) 

= -g&'[A li ,A v ]+g[F llv ,A v ]. (3.8) 

and G(x,y) = — (x\d~ 2 \y). As we have written it, the term in ( p.7| ) involving ghosts cannot 
immediately be replaced by derivatives on JF. For this, we proceed as follows. Integrating 
out the ghost fields in T we get Faddeev-Popov determinant det(-d-D) = e tr ln (~ 9 D ) ; which 
is equivalent to a term in the action — ztr ln(— d ■ D). This leads to a ghost-current of the 
form 



5A« 



oo 

71=1 



Y^(-g) n+1 J Tr pTM^O ■ ■ •^( 2 / n )] a *-^G(a; j3 / 1 ) ■ • • d Vn G(y n ^, y n )d„G(y n , x) , 

(3.9) 



where we have introduced the adjoint representation of T a 's by (T^ d j) bc = —f abc . (C^ t in 
this notation reads Cf^t — ~9diicFadj c ^i-) 

The trace in the above equation may be written in terms of the generators in fundamental 
representation using f abc = — 2Tr T a [T b ,T C ], 

Tr (T ai T ai ■ ■■T an ) adj = 2Tr T a [T\ [T a \ [T a \ [. . . , [T a "~\ [T an , T b ]] . . .] . (3.10) 

Therefore we can finally write the ghost part of the current (as a matrix in the fundamental 
representation) as 



J,(A) gh = j;{A) gh T a (3.11) 

oo 

iJ2(-9) n+1 D^...» n (x,yn ■ ■ -,yn)[T b , [A^( yi ), [..., [A^iy^), [A^(y n ),T b ]) . . .]. 

71=1 



where 



Dp ftl ...f ln (x,yi, ...,y n ) = d^Gfayx) ■ • •d (tn G(y n .. 1 ,y n )d (t G(y n ,x). 

Then ( |3.7|) finally becomes 
.5T 



(3.12) 



A=iG4- 



(3.13) 



For connected part C[a] = — zln.F[a], 



5C 



n— = d 2 a a u + J a AA) + hJ a (A) 
6 a? 



gh 



A=-hG s £+ihG£- 



(3.14) 



This is the basic equation for the scattering amplitudes. The n-expansion of this equation 
leads to recursion rules upto arbitrary order. We have explicitly displayed the factors of h in 
( |3.14j ). Recall that the interaction part of the action carries 1/h and that each propagator 
G carries K. The ghost terms, arising from the determinant, have h°. C[a] in ( |3.14j ) has an 
expansion of the form 



C[a] = -iln ^[o] = ^C (0) + + hC^> + 



(3.15) 



Starting from (|3.14 ) and using the above expansion for C[a], one can, in principle, 
systematically derive recursion relations to any desired order in h. 

Let us begin with a tree- level recursion formula for . This is simply given by 

6C® 



(3.16) 



where 



A(0)a _ 



G 



5C<® 



(3.17) 



i.e., 



d*Af = d*a» + JMf)- 



(3.18) 



This is just the equation of motion as it should be. is essentially the same object as 
the current that Berends and Giele used to derive tree-level recursion relation for gluon 
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scattering processes; if one expands in powers of a M , the coefficient function of each term 
gives one-gluon off-shell current of Ref. [Q] when multiplied by polarization vectors of on- 
shell external gluons. More generally, if one considers a^'s off shell, it gives the generalized 
current with off-shell gluons which has been used for qq —>■ qqgg ■ ■ ■ g process and some one- 
loop calculations |lT],|12j. Also, notice that from the form of current J M in (|3.8|) , it is obvious 



that color factors are factorized from coefficient functions and we can write 

oo „ 

4 0) = E / CgL^fo Vi, • • • ,Vn)a^(yi) ■ ■■a^{y n ) . (3.19) 

n=l J 

The coefficient functions C^'s do not carry color indices. We also have Cjj ) / } l (x,yi) = 
S{x - yi )5^. Using (£T3) in ( fOg) , 



d 2 E / C®(x, Y)a(Y) - d 2 a, = J,(A) 

= -E/ [WCf{x,Y)Cf\x,Z) - 2d u C^\x,Z)C^\x,Y) + Ci \x,Y) d, C< \x,Z) 
+C^(Y)d u C' /(0 \x,Z) - Cf\x 1 Z)& / C^\x,Y)]a{X)a{Y) 
-E / [2C?>(*,A:)<7< ^^ 

-C^(x,X)Ci°\x,Y)C u{0) (x,Z)}a(X)a(Y)a(Z) , (3.20) 

where X, Y and Z are collective indices and a(Y) stands for ]li <Vi(^«) e ^c. If a M is restricted 
to be on-shell, the terms containing d u C^ in the third line vanish because of current con- 
servation. We can transform this equation to momentum space by writing 

C^ 1 ..., n (x,y 1 ,...,y n ) = /(27r) 4 (* + Eft)^..^(Pi,---^)e*"e <Eft *, (3-21) 
where momentum conservation has been taken into account. Then ( |3.20 ) becomes 



n-1 



Pl n C^\l, ...,n)= VW,m, P m+ x,n)Ci°\l, m)Cf\m + 1, . . . , n) 



m=l 

n— 2 n— 1 



+ E E V i r a Ci°\l,...,m)C^(m + l,...,k)C^(k + l,... 

m=l k=m+l 



(3.22) 



where Py — Pi + ■ ■ ■ + Pj and V pup , V pupa are color-ordered vertices, 



vf vp ( P , q) = -g{g vp {v - qT + g p "(2q u + f) - gTW + <f)}, 

vr pa = -g(2g"gr - g pv g pa - g^g va ) • (3.23) 



This equation is the recursion relation for currents with off-shell gluons derived in B, [O 



In deriving higher order recursion relations for S-matrix elements, we must correct (|3.14) 



by including the appropriate renormalization constants. As usual, we interpret fields and 
coupling constants as renormalized ones and assume that renormalization counterterms are 
included in C int of (|3.4j) . Explicitly, 



C mt = ~6Z 3 Tr (d,A v - d u A^f - 2gZ.Tr d^A v [A^ M - ^ZfZ^Tr [A^ A v f 

-25Z 3 Tr c(-d 2 )c - 2gZ 1 Z 3 /Z 3 Tr d„c{A^ c] , (3.24) 

where y/Z^A^, \fz~ 3 c and g Z\Z~ 3 3 ^ 2 are bare quantities with 5Z 3 = Z 3 — l,5Z 3 — Z 3 — 1. 
Correspondingly, becomes 

Jp(A) = —gd v [Afj,, A v ] + g[F^ U) A v ] 

~5Z 3 d v {d„A v - d u A„) - g{Z 1 - l)(d v [A», M - [F^, A v \) 

+g 2 Z. {Z x /Z 3 -1)[[A^,A U }, A v ] . (3.25) 

The extra counterterms contribute to one-loop or higher orders and will cancel ultraviolet 
divergences from loop integrals. In addition, to account for the freedom of arbitrary finite 
renormalization we include a finite wavefunction renormalization constant z 3 in the external 
field, i.e., a^^a^ = a^/^/z^. 

Now we are ready to discuss one-loop recursion relations. Keeping terms relevant upto 
one-loop order, 



An = —G— G— ^z 3 G— — + ■ 

4 0) +4 1) + --- ( 3 - 26 ) 



The current can similarly be expanded as J^{A) = + + • • ■, with 
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4 0) 



-gd v [Af,AW]+g[F» v {A {0) ),^>], 
-gd v {[Af,A^\ + [A^\A^\)+g[F, v {A m \A^\ 

sz 3 . 



g[D„AU-D v AW,A^ 



(0)N 



'-d»{d,Af -d u Af) 



(3.27) 



where D^ = d^ + g[A@\ }. In ( gig ) we need 



J„(i4 + iG^)i = j M (i4) - gSTW^ 4 0) ] + sK G ^r/ - ^ 0) ] 



(3.28) 



Upto the one-loop order, the functional derivative term in the ghost current has no contri- 
bution, 



(3.29) 



J,{A + iG—) gh l = J^) gh + ■ ■ ■ . 
Collecting all these, we get an equation for A^\ 

d 2 AU = jU + UA ( % h 

-gd u [tG^-, A^] + g[td^G^- - id u G-^—, A®] 

~9 2 [iVG^AflA^] + [[A$),i G A-],A^ 

where A^ is the solution of ( |3.18|) . Notice that, in this case, iGjj^ — terms contract color 
indices when acting on A^ and so the color decomposition does not occur as in tree-level 
case. However, it is possible to write A^> as a sum of color-factorized amplitudes and the 
proof has been given using string-theory argument [13] and color flow diagrams [|16|| . Here 
we give a simple proof based on (I 



(3.30) 



First we study the action of iG-£- on A<£\ From (|3TTP| ), 

^ oo n „ 

[Gt-7,4 0) ] = 53 53 / c il\..M m .^u m+1 ..M n (^,yi, ■ ■ ■ ,yrn-i,y,y m +i, ■ ■ ■ ,y n )G(x,y) 

n=l m=l 

x [T a , a^(yi) • • • a I/m _ 1 (y m ^i)T a a Um+1 (y m+1 ) • • • a Vn (y n )\ . (3.31) 
With the help of Fierz identity for SU (N) 
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(T a XT% = - l - (SyTt X - , (3.32) 



it becomes 

n=l m=0 



[Gj^ } 4 0) ] = EE / 1, • • • , m; m + 1 . . . , n)o(l) • • • a(m)Tr [a(m + 1) • • • o(n)] , 

(3.33) 



where 



~ If 
C Ufl {x,l,...,m;m + l...,n) = -- / G(x, y)[(7g u 



~Gw m ...v n ^l...v m { X -i • • • ) 2/n) 2/) 2/1) • • • ) 2/m)] • (3.34) 



(Here the 1/iV-term in ( p. 32 ) does not contribute because of the commutator; notice also 



that the summation in (|3.33 ) starts with m = 0.) Thus differentiation of A^ a produces 



terms with trace over substrings of a^'s. Also, from ( |3.11 ), we see that the ghost current 
has the same structure, 

oo 

n=l 
n 

x £ £{ T ^ (0) (*i) • • • 4 (0) fe)T^( %1 ) ■ • • A^(i n ) 
k=o {n} 

-A® fa) ■ ■■A^(i k )T h A^(i k+1 ) ■ ■■A^( ln )T b }, (3.35) 

where is the sum over all permutations of {l,...,n} such that i\ < ... < i k and 

i k+ i > ... > i n . Using ( |3.32|) , we get 



i 00 n 



^% = -«EE E ^ +1 (-i)^^i>---^n) 

^ n=l fc=0 o-e5 n;fc 

xA(°)(l)---AW(A;)Tr [#(fc + l)---A (0) (n)], (3.36) 



where 5 n; fe is the set of all permutations of {1,2,..., n) that preserves the ordering of {a} = 
{1, . . . , k} and the cyclic ordering of {/3} = {n, n—1, . . . , k+1}, while allowing for all possible 
relative orderings between elements in the two sets (For example, (1, n, 2, . . . , k, n — 1, . . .) 
is in S n . k but (2, 1, . . . , k, n, n — 1, . . .) is not.); is given by 
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D^x, 1, . . . , n) = D^x, 1, . . . , n) - (-l) n D^x, n, n - 1, . . . , 1). (3.37) 
Since (|3.30|) is at most linear in A^\ it is clear that can be written as 

oo n „ 

A< £H X ) = £ £ y CSi^) l,...,m;m + l,...,n)a(l)---a(m)Tr [a(m + 1) • • • a{n)] . 

(3.38) 



n=l m=0 ' 



Then in flOoD we can separately equate terms with the same trace structure and 's 



with different m's do not mix with each other. The functions we need are then = 
(x, 1, . . . , n) corresponding to the term in (|3.38|) with no trace. These obey the recursion 
rule obtained by substituting = J2™=i I C^(x, 1, . . . , n)a(l) ■ ■ • a(n) in (|3.30| ). The other 
amplitudes can be obtained from CjQ. 

The coefficient functions corresponding to different trace structures have simple relations 
among them (which allow us to construct the amplitudes with subtraces from C^a). We 
have 0, 



CSJ i (x,l,... ) m;m + l ) ...,n) = (-l)"- ra £ C^{x,a u . . . ,a n ) , (3.39) 



For this, we show that the right hand side of ( |3.39|) satisfies the same equation as that for 



the left hand side. This is most easily seen for the ghost current which we shall consider first. 
Obviously, an equation like ( |3.39|) holds for D^s (with the identification of the summation 
indices k in (|3~36p and m in (|3~39|) ), if Aj®{x) = JCj®(x,X)a(X) in (|3~36|) is replaced 



by its lowest order term a^. For the terms with more than one a^'s from one A^°\ there 
are potential discrepancies between both sides. This is because in ( |3.39| ) the sumation is 



over all o G S n - m while, in Q3.36| ), we sum only over a G S n -k (k < m) which does not 



include such permutations that mix indices from both {a} set and {/3} set within one tree 
structure. However, the extra terms in fl3.39| ) cancel out thanks to the symmetry properties 
of color-ordered vertices, 

ypU P * + ygpVO + yJUTUp = Q ? (3 ^ 
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which, when applied to ( |3.22|) together with induction, leads to identities for 1 < m < n, 

cfK...,og = o, (3.4i) 

where S' n . m is the same as S n - m defined above except that it does not include the cyclic 
permutations of {f3} set. [This equation can be considered as a generalization of the so- 
called "dual Ward identity" for tree amplitudes which corresponds to m = 1 case |§, [[17]] •] 
Thus (|3.39|) holds for ghost current. 



We shall now show that the relation ( |3.39|) connecting amplitudes with different trace 
structures holds for the nonghost terms in ( |3.38| ) as well. Towards this, consider the terms 
with differentiation G-^ in (|3.30|) . Diagramatically, those terms correspond to one-loop 



gluon diagrams with the external leg x directly connected to the loop. In our setting, we 
can proceed as follows. Applying (|3.22| ) repeatedly to Cfi'(x, 1, . . . , m, y, m + 1, . . . , n) we 
can identify the internal line for each term from ( |3.22| ). It connects x and y and the other 
external legs are ordered in clockwise direction in such a way that 1, . . . , m-th legs are below 
the line while m + 1, . . . , n-th legs are above the line (Fig. 1). Then we draw a line connecting 
x and y which enclose m + 1, ... , n-th legs so that it represents that the legs inside the loop 
are traced. Then, it is easy to see which terms generate which color structures. Obviously, 
given such a diagram, we get the same subtrace structure from diagrams made by altering 
the relative order of trees belonging to different sets while keeping the order of {1, . . . , m} 
and the cyclic order of {m + 1, . . . , n} (Fig. 2). Moreover, they are trivially related to terms 
which contribute to C^L i.e., we can simply pull out the trees inside the loop to outside with 
minus sign, using symmetry property of vertices ( [3.40| ). Notice that during this procedure 
the order of {m + 1, . . . , n} is also reversed. Thus, essentially we have the same situation as 
in the case of ghost current and the relation of the type ( |3.39| ) holds in this case as well. 

Now it remains to consider the term from which contains C}}> 's with less number of 
legs. It corresponds to the gluon-loop diagram with the external leg x attached to a tree. 
This case, however, is not much different from the previous one and we can argue in the 
same way with the help of ( |3.41|) if needed. Finally, there are counterterm contributions 
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for CW(s, 1, . . . , n) case in contrast to C^} (x, 1, . . . , to; to + 1, . . . , n) (to < n). But again 
the identity (|3.41|) guarantees that those terms cancel out when the summation in (|3.39|) is 
done. This completes the proof that both the left and the right hand sides of (|3.39|) satisfy 



the same recursion relation. Since (|3.39|) trivially holds for n = 3, they are indeed equal to 
each other. 

Beyond one-loop, it is clear from ( |3.14| ) that A^' will in general have terms with k 



subtraces, 

oo - 

A ^\ x ) = J2 I C m! m2 ...m k ^, 1 ,---,m 1 -,m 1 + l,...,m 2 ]...- 1 m k + l...n) 

n=1 {mi} 

xa(l) • • • a(mi)Tr [a(mi + 1) • • • a(m2)] • • • Tr [a{rrik + 1) • • ■ a(n)} , 

(3.42) 

and each C^ ma< mkfM with different k will satisfy its own equation. It might be possible to 
find simple relations between them as in one-loop case. 
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FIGURES 

FIG. 1. A typical term in C^\x, 1, . . . , m, y, m + 1, . . . , n). y is connected to x and the loop 
encloses legs m + 1, . . . , n which are under a trace. 

FIG. 2. (a) A diagram obtained by changing the relative order of traced legs and the others 
from Fig. 1. It contributes to Cm^(x, 1, . . . , m; m + 1, . . . , n); (b) A diagram obtained by pulling 
the legs out of the loop in Fig. 2(a). It contributes to Cn1[x, 1, . . . , n). 
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Fig. 2 (b) 



